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Abstract 

The transverse momentum (px) dependence of the nuclear modification factor /?aa and the centrality 
dependence of the average transverse momentum (px) for inclusive ]/\j/ have been measured with 
ALICE for Pb-Pb collisions at = 2.76 TeV in the e+e^ decay channel at mid-rapidity (|y| < 0.8). 
The (px) is significantly smaller than the one observed for pp collisions at the same centre-of-mass 
energy. Consistently, an increase of /?aa is observed towards low px. These observations might be 
indicative of a sizable contribution of charm quark coalescence to the J/xj/ production. Additionally, 
the fraction of non-prompt J/y/ from beauty hadron decays, /b, has been determined in the region 
1.5 < Px < 10 GeV/c in three centrality intervals. No significant centrality dependence of /b is 
observed. Finally, the 1 ?aa of non-prompt J/i/ is discussed and compared with model predictions. 
The nuclear modification in the region 4.5 < px < 10 GeV /c is found to be stronger than predicted 
by most models. 
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1 Introduction 


Heavy-ion collisions at high energies allow the study of strongly interacting matter under extreme condi¬ 
tions. Calculations based on Quantum-Chromo-Dynamics (QCD) on the lattice indicate that the hot and 
dense medium created in these collisions behaves like a strongly coupled Quark-Gluon Plasma (QGP) 
M. Heavy quarks are an important probe for the properties of this state of matter, since they are 
produced via hard partonic collisions at a very early stage and thus experience the complete evolution of 
the system. Quarkonium states, i.e. bound states of a heavy quark and anti-quark such as the J/i/r meson 
(cc state) are of particular interest. It was predicted that the 3/y formation is suppressed in a QGP due to 
the screening of the cc potential in the presence of free colour charges llsl]. Experimentally, a suppression 
of the inclusive J/i/r yield in heavy-ion collisions relative to the corresponding yield in pp, scaled by 
the number of binary nucleon-nucleon collisions, has been observed at the Super Proton Synchrotron 
(SPS) M and the Relativistic Heavy Ion Collider (RHIC) fllli. The level of suppression was found 
to be similar at SPS and RHIC, despite the significantly different collision energy. More recently, the 
nuclear modification of J/i/r was also measured for Pb-Pb collisions at the EHC Ill-l3l- While at high 
transverse momentum (pj > 4 GeV/c) the suppression factor is at the same level as the one observed 
at RHIC in the low pj region, a significant reduction of the suppression is measured towards lower pj. 
This has been interpreted as the effect of an additional contribution to J/i/r production at low pT, due 
to the combination of correlated or uncorrelated c and c quarks il4l IlSll. This contribution becomes 
sizable at EHC energies, since the number of cc pairs is much higher than at lower energies. Assuming 
that a deconfined phase is produced and that all the J/i/r are dissociated, this process happens at the 
chemical freeze-out stage of the fireball evolufion. This is fhe approach followed wifhin fhe sfafisfical 
hadronizafion models described in Refs. iSEl. Alfernafively, J/i/r could be generated via coalescence 
fhroughouf fhe full evolufion of fhe QGP phase, if fheir survival probabilify in fhis environmenf is large 
enough. This scenario has been implemenfed in several parfonic fransporf models |1^ 19]. It was found 
that both approaches can provide a description of the measured nuclear modification factors il2ll and of 
the elliptic flow of inclusive J/xj/ 1201] . 


The production of open beaufy hadrons is expected fo be sensifive fo fhe densify of fhe medium created 
in heavy-ion collisions due fo fhe energy loss experienced by fhe parenf parfon (a beaufy quark) which 
hadronizes info fhe beaufy hadron. This energy loss is expected fo occur via medium-induced gluon ra¬ 
diation i2lL l22ll and elasfic collisional energy loss processes 123142511 and if depends on fhe QCD Casimir 
coupling faclor of fhe parfon (larger for gluons fhan for quarks) and on fhe parfon mass ll26l429ll. Ofher 
mechanisms, such as in-medium hadron formafion and dissociation, can be envisaged as particularly 
relevanf for heavy-flavour hadrons due fo fheir small formafion times i30l432n . 


Inclusive 3/y producfion is fhe sum of several confribufions. In addition fo fhe direcfly produced J/i/r, 
fhe decays of heavier charmonium sfafes, such as fhe Xc and i/r(2S), also confribufe fo fhe inclusive 
3/y yield. These fwo sources (direcf and charmonium decays) are defined as prompf 1/ 1 //, where fhe 
contribution from charmonium decays is about 35% as measured in pp collisions ll33n . Since heavier 
charmonia are less strongly bound than the 3/y they should be more easily dissolved in a deconfined 
medium 041] . The 3/\l/ suppression measured af fhe SPS is indeed compatible wifh fhe assumpfion fhaf 
only fhe excited sfafes are dissolved and nof fhe direcfly produced J/i/r flQj. On fop of fhe prompf J/i/r 
producfion, fhere is an additional non-prompf confribufion fo fhe inclusive J/i/r af high cenfre-of-mass 
energies, coming from fhe decay of beaufy hadrons. Since fhese decays proceed via weak inferacfions, 
fhe resulfing 3/y will originafe from a decay verfex thaf is displaced from fhe main inferacfion vertex. 
Their measuremenf provides a direcf determination of fhe nuclear modificafion of beaufy hadrons. By 
subfracfing fhe non-prompf confribufion from fhe inclusive 1/y/ yield one can also provide an unbiased 
information on medium modification of prompf charmonia. The non-prompf 3/y confribufion at mid¬ 
rapidity has already been measured in pp collisions at ^/s = 1 TeV by ATEAS ll35n . CMS Bal and 


AEICE OTil . Eor Pb-Pb collisions at 


' = 2.76 TeV CMS has also published prompt and non-prompt 
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J/i/r production results at mid-rapidity for pT >6.5 GeV/c il3h . 


In this paper we present a differential measurement of the inclusive 1/ ip production at mid-rapidity in 
“ = 2.76 TeV. The pj dependence of the nuclear modification factor and the 


Pb-Pb collisions at 


centrality dependence of the average transverse momentum of J/ i/r have been obtained, extending the set 
of results presented in lll2l] . A measurement of the prompt and non-prompt contributions to the inclusive 
J/i/r production is also presented. The nuclear modification factor of non-prompt J/ip is determined 
down to pt = 1-5 GeV /c and compared to model predictions. 


2 Data Analysis 


A detailed description of the ALICE detector can be found in ll38ll . Eor the analysis presented here the 
detectors of the central barrel have been used, in particular the Inner Tracking System (ITS) and the 
Time Projection Chamber (TPC). These detectors are located inside a large solenoidal magnet with a 
field strength of 0.5 T. They allow the measurement of J/yr mesons via the dielectron decay channel 
in the central rapidity region down to zero px. The ITS 091] consists of six layers of silicon detectors 
surrounding the beam pipe at radial positions between 3.9 cm and 43.0 cm. Its two innermost layers 
are composed of Silicon Pixel Detectors (SPD), which provide th e sp atial resolution to separate on a 
statistical basis the non-prompt J/y/. The active volume of the TPC IITOll covers the range along the beam 
direction —250 < z < 250 cm relative to the Interaction Point (IP) and extends in radial direction from 
85 cm to 247 cm. It is the main tracking device in the central barrel and is in addition used for particle 
identification via the measurement of the specific ionization (dE/dr) in the detector gas. 


Triggering and event characterization is performed via forward detectors, the VO 114ih and two Zero 
Degree Calorimeters (ZDC) 14211 . The VO detectors consist of two scintillator arrays positioned at z = 
—90 cm and z = -1-340 cm and cover the pseudo-rapidity ranges —3.7 < rj < —1.7 and 2.8 < rj < 5.1. 
The ZDCs, each one consisting of two quartz fiber sampling calorimeters, are placed at a distance of 
114 m relative to the IP in both directions along the beam axis and are used to detect spectator nucleons. 

The results presented in this article are based on data samples collected during the Pb-Pb data taking 
periods of the LHC in the years 2010 and 2011. In the case of the 2011 data sample the Minimum Bias 
(MB) Level-0 (LO) trigger condition was defined by the coincidence of signals in both VO detectors along 
with a valid bunch crossing trigger. Lor the 2010 data sample, in addition, the detection of at least two 
hits in the ITS was required. Both MB trigger definitions lead to trigger efficiencies larger than 95% for 
inelastic Pb-Pb collisions. Electromagnetic interactions were rejected by the Level-1 (LI) trigger, which 
required a minimum energy deposition in the ZDC by spectator neutrons. The beam-induced background 
was further reduced during the offline analysis by selecting events according to the relative timing of 
signals in VO and ZDC. The offline centrality selection is done using the sum of the two VO signal 
amplitudes. By fitting the corresponding distribution with the results of Glauber model simulations, the 
average number of participants (Apart) and the average nuclear overlap function (T^a) = (Acoii) /for 
a given centrality class can be determined as described in J^ . Here, (Vcoii) is the average number of 
binary nucleon-nucleon collisions and the inelastic nucleon-nucleon cross section. The numerical 
values for (Apart), (Acoii), and (Taa) are tabulated in 112]. 

The 2010 data sample consists of 1.5 x 10^ events, taken with the corresponding MB trigger. The 2011 
event sample was enriched with central and semi-central Pb-Pb collisions by using thresholds on the 
VO multiplicity at the LO trigger. Lrom the latter data set we analyzed 1.9 x 10^ central (0-10% of the 
centrality distribution) and 1.7 x 10^ semi-central (10-50%) events. The summed 2010 and 2011 data 


samples correspond to an integrated luminosity of = 26.4 ± 0.3(stat.)^j y(syst.) ph * lll2ll . 
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Fig. 1 ; The invariant mass distributions of inclusive J/y/ at mid-rapidity (|y| < 0.8) for Pb-Pb collisions (0-40% 
most central) at = 2.76 TeV. The left panels show the interval 0 < pj < 2.5 GeV/c and the right ones 
2.5 < pt < 6 GeV/c. The upper panels display the opposite sign distributions together with the result of the 
mixed event procedure. In the lower panels the background subtracted distributions are shown and compared to 
the simulated line shape. Also, the signal-to-background ratio S/B and the significance of the signal are given. 


2.1 Inclusive J/i/r 

J/Y candidates are reconstructed by combining opposite-sign (OS) pairs of electron/positron candidates 
and calculating their invariant mass mge- These candidates are selected from tracks reconstructed in the 


ITS and the TPC by employing the set of quality criteria described in il2L 14411 . In order to reject the 
background from photon conversions in the detector material, tracks are required to have a hit in one 
of the SPD layers. In addition, at least 70 out of a maximum of 159 space points reconstructed in the 
TPC must be assigned to a given track, which also needs to fulfill a quality criterion of the track fit 
(X^/ndf < 4). The tracks are required to be in the range |t]| < 0.8, where the tracking and particle 
identification performance of the TPC is optimal, and to have pj > 0.85 GeV/c to improve the signal- 
to-background ratio in the 1/ y/ mass region. 

Electron candidates are selected by requiring that the dE/dx measurement in the TPC lies within a band 
[—la,-|-3a] around the momentum-dependent parameterization of the expected signal, where a is the 
phase space dependent dE/dx resolution (details can be found in 14511 ). The selection is asymmetric in 
order to minimize the contribution from pions. To further suppress the hadron contamination, tracks that 
are compatible within zh4a with the proton expectation are rejected. A side effect of this cut is that tracks 
below pt = 1 GeV/c are effectively removed. 
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Measurement of the inclusive J/^f yield 

The J/ y/ signal counts Ajare obtained from the number of entries in the background subtracted invari¬ 
ant mass distributions in the range 2.92 < ntge < 3.16 GeV/c^. The uncorrelated background is evaluated 
with a mixed event (ME) technique. In order to achieve a good description of the background only elec¬ 
trons and positrons from events with similar properties in terms of centrality, primary vertex position, 
and event plane angle are combined. The ME distributions are scaled to the same event (SE) distribu¬ 
tions in the mass ranges 1.5 < < 2.5 GeV/c^ and 3.2 < wige < 4.2 GeV/c^, so that the J/yr signal 

region is excluded. The normalization area contains the y/^(2S) signal, but its contribution is negligible 
and can therefore be safely ignored. Also, contributions from the tail of the J/yr signal shape to this mass 
interval are below the percent level and will thus not significantly affect the normalization. Eigure [T] 
shows a comparison of the SE and ME invariant mass distributions for the 0-40% most central Pb-Pb 
collisions for electron-positron pairs at mid-rapidity (|y| < 0.8) in two px intervals: 0 - 2.5 GeV/c and 
2.5 - 6 GeV/c. The agreement between the SE and ME distributions outside the signal region is very 
good and allows signal extraction with significances larger fhan eighf. 


The J/yr yield per MB even! in a given px interval, Tj/y,, is obfained as 


Yj/y{Pt) = 


^j/v/(pt) 


BRgg Agvts (Axe)(px) ■ 

Here BRgg is fhe branching rafio for fhe decay J/y/ —)• e+e^, Agvts the number of evenfs, and (A x e) fhe 
phase space dependenf producf of accepfance A and reconsfrucfion efficiency £. The laffer is calculafed 
from Monfe Carlo (MC) simulafions as fhe rafio befween fhe number of reconsfrucfed and generafed MC 
1/yr, which are assumed fo be unpolarized. In pp collisions af ^/s = 1 TeV fhe J/y/ polarizafion has 
been measured and was found fo be compafible wifh zero af mid-rapidify (px >10 GeV/c) and forward 
rapidity (px > 2 GeV/c) ll46l-t48n. In heavy-ion collisions no measuremenf exisfs, buf J/y/ mesons 
produced from fhe recombination of charm quarks in fhe medium are expected fo be unpolarizecOl The 
MC evenfs used for fhe calculafion of (A x e) are consfrucfed by adding fo background evenfs, generafed 
wifh fhe HIJING model lH^, J/y/ mesons decaying info e+e^ pairs, whose phase space disfribufion 
is obfained from exfrapolafions of ofher measuremenfs ll^ . faking info accounf shadowing effecfs as 
parameterized in EKS98 15111 . The dielecfron decay is simulafed wifh fhe EvfGen 115211 package, using 
fhe PHOTOS model 15311 fo describe fhe influence of final sfafe radiafion. This choice, fogefher wifh fhe 
simulafion of bremssfrahlung in fhe defecfor material, is mandatory for a proper descripfion of fhe low 
mass fail in fhe measured J/y/ mass disfribufion and ensures fhaf fhe fracfion of fhe signal oufside of 
fhe /Wee infegrafion window is properly accounfed for in fhe correcfion (Axe). The propagafion of fhe 
simulafed particles is done by GEANT3 15411 and a full simulation of fhe defecfor response is performed. 
The same reconsfrucfion procedure and cufs are applied to MC evenfs and fo real dafa. The qualify of fhe 
simulafion is illusfrafed by fhe good agreemenf of fhe background-subfracfed invarianf mass disfribufions 
wifh fhe MC simulafion of fhe J/y/ signal shape, after normalizing if to fhe same infegral as fhe measured 
signal (see Pig.Ull. 


The analysis has been performed in fwo slighfly differenl cenfralify infervals (0-40% and 0-50%), where 
fhe larger one is used for fhe exfraclion of non-prompf J/y/ which requires a higher sfafisfics fhan fhe 
inclusive measuremenf. Also, fhe px infervals have been optimized for fhe differenl analyses. If was 
checked fhaf fhe resulls for inclusive J/y/ obfained wifh fhe fwo cenfralify binnings are in good agree- 
menl. 


Determination of the pp reference for Rj^j^ 

Prom fhe corrected J/y/ yield Tj/,^(px) the nuclear modification factor Raa(pt) is calculated as 


Raa(Pt) 


(T’aa) ■ 


* The impact of the polarization on the acceptance was studied for extreme polarization scenarios in m. 


(2) 
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Fig. 2: The average transverse momentum (pj) of e+e^ pairs, measured for the pj range 0-10 GeV/c, as a func¬ 
tion of the invariant mass mge in centrality selected Pb-Pb collisions at = 2.76 TeV. The shown uncertainties 
are statistical only. The background (pj) distributions and the total fit results are also shown superimposed to the 
data points. 


Since no differential J/i/r measurement at mid-rapidity at low pT is available for pp collisions at ^/s = 
2.76 TeV 1551] . the reference needed for the construction of /?aa is based on an interpolation of the mid¬ 
rapidity measurements by PHENIX at ^/s = 0.2 TeV |56], CDE at ^/s = 1.96 TeV |57], and AEICE at 
^/s = 7 TeV 1155n . The interpolated pj distribution is obtained by fitting the following parameterization 
to the available data sets Isoll 

1 d% a PJ 


= c ■ 


da/dy dzxdy (1 + a^Zj) 

Here, zx is defined as pt/(pt)> ^ = r(3/2) T{n — 3/2)/Y{n — 1), and c = 2{n — 1) a, where n is 
fhe only free fil paramefer. The value for (px) (calculated in fhe px range 0-10 GeV/c) af y/s = 
2.76 TeV, which is needed fo franslafe fhis parameferizafion info da/dpx, is determined by inferpo- 
lafing befween fhe exisfing (px) measuremenfs for pp and pp collisions 113-57]. This interpolation is 
done using various functional forms for fhe ^/s dependence fo defermine fhe systematic uncerfainfy. 
Eor fhe absolute normalization of fhe para mefrized specfrum, fhe same interpolated value Ao/Ay = 
4.25±0.28(sfal.) ±0.43(sysf.) pb as in |12] is used. 


The main sources of sysfemafic uncerfainfies for fhe px dependenf 7 ?aa of inclusive J/y/ are fhe signal 
reconsfrucfion procedure, fhe MC inpuf kinemafics, fhe uncerfainfies on fhe inferpolafed pp reference 
and on fhe nuclear overlap funcfion. The corresponding values are summarized in Tab. [T] White fhe firsl 
fwo componenfs are uncorrelafed befween fhe px intervals (type II), the uncertainty due to the nuclear 
overlap function is fully correlated (type I). The pp reference on the other hand introduces both uncorre¬ 
lated and correlated contributions. To determine the uncertainty related to the signal reconstruction, the 
normalization range of the ME background and the size and positions of the mgg bins have been varied. 
All track and electron selection criteria, such as the electron inclusion cut and the SPD hit requirement, 
have been relaxed and/or tightened in order to test the stability of the result, as was performed in 11211 . 
The value of the systematic uncertainty is determined as the standard deviation of the distribution of all 
results obtained with the listed variations. The evaluation of the uncertainties associated with the MC 
input kinematics is also described in tfl2l] . white the uncertainty of the pp reference is estimated from the 
differences between the cross-section values obtained with the fitting procedure based on Eq. ([3]) and the 
measured values used for the fit at the various energies. 


Determination of (px) and {p\) 

Since the collected Pb-Pb statistics would allow the extraction of the J/i/r yield in a few px intervals only, 
the average transverse momentum (px) is determined by a fit to the distribution of the (px) of e+e^ pairs 
as a function of rngg. When building such a distribution, the individual e+e^ pairs are weighted by the 
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Pt range (GeV /c) 




0^0% 


0-50% 



Source 

0-2.5 

2.5-6 

0-1.5 

1.5-4.5 

4.5 - 10 

Type 

Signal reconshuction 

10.9 

11.0 

22.6 

8.5 

10.0 

II 

MC input kinematics 

3.2 

5.9 

1.5 

A.l 

5.3 

II 

Nuclear overlap function (Taa) 

3.2 

3.2 

3.8 

3.8 

3.8 

I 

pp reference 

12.1 

12.1 

12.1 

12.1 

12.1 

I 

4.5 

4.9 

5.0 

4.3 

10.1 

II 

Total 

12.6 

12.6 

\1.1 

12.7 

12.7 

I 

12.3 

13.5 

23.2 

10.6 

15.2 

II 


Table 1: The correlated (type I) and uncorrelated (type II) systematic uncertainties (in percent) on the measure¬ 
ment of the nuclear modification factor /?aa of inclusive J/y/ for several pj intervals in Pb-Pb collisions (0-40% 
and 0-50% most central) at = 2.76 TeV. 

inverse of their acceptance times efficiency (A x assuming that they come from the decay of a J/i/r. 
The resulting (pj) distributions are fitted by the expression 

(PT)meas = \ [“^Ke) (Pt) J/y. + (PT)Bkg] ■ (4) 

O \ fnQQ j -\- D [JtLqq j 

Both factors S and B depend on mge and correspond to the distribution of the J/i/r signal and of the 
background. For S the same background subtracted signal distribution S{mes) is used as for the extraction 
of the yield (see lower panels of Fig. [T]), while the background B is generated from the ME sample, as 
Bim^e) = c^B^^{mss). The normalization factor is determined by fitting to the corresponding 

/Uge distribution of e+e^ pairs in the regions 1.5 < rngg < 2.5 GeV/c^ and 3.2 < rngg < 4.2 GeV/c^, thus 
excluding the signal region. For the sum S(mee) -|-B(/7iee) in the denominator of Eq. Q, the measured 
OS pair rngg distribution is used. The (pT)Bkg> defined as the (pj) of the combinatorial background pairs, 
is also calculated from the ME sample. This analysis is performed in three different centrality intervals: 
0-10%, 10^0%, and 40-90%. Figure |2] shows the measured (pj) of the e+e^ pairs in the pj range 
0-10 GeV/c together with the results of the fit procedure. In addition, with an equivalent method, the 
mean square transverse momentum {p\) is also calculated for the same centrality intervals. 

The systematic uncertainties of the (pt) measurement for inclusive J/ i/r are mainly determined by the 
signal extraction, the stability of track and electron selection criteria and the fit procedure (see Tab. O. 
While the first two components are not correlated between the different centrality intervals (type II), the 
systematic uncertainty intrinsic to the fit procedure can affect the data points in a correlated way (type I). 
The uncertainties related to the signal extraction have been evaluated by varying the normalization range 
of the ME background, the size and positions of the rngg bins, the fit region and by using in addition 
to the ME background a linear function for the background description. In addition, the stability of the 
fit procedure was tested by modifying the approach, e.g. by using in the fit the direct sum of S and B, 
instead of the OS pair distribution, or by using a fit function for B, instead of ME. It was also verified 
by applying the above described method to MC events, which were constructed by combining signal 
with background events with a realistic S/B ratio. It turned out that the procedure allows to recover the 
(pt) of the simulated J/y/ mesons within a 2% difference. This value is assumed as correlated (type I) 
uncertainty. Finally, the uncertainty in the signal-to-background ratio is propagated into the statistical 
uncertainty of the (pt)j/v/- 

2.2 Non-prompt J/i/r 

The candidate selection for the non-prompt J/i/r analysis includes, in addition to the previously described 
criteria, the condition that at least one of the two decay tracks has a hit in the first SPD layer, in order to 
enhance the resolution of secondary vertices. 
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Source 

0-10% 

Cenfralify 

10^0% 

40-90% 

type 

Signal exfracfion 

2.0 

1.6 

4.5 

II 

Track selecfion 

3.1 

3.9 

13.7 

II 

Sfabilify of fhe fif procedure (from MC) 

2.0 

2.0 

2.0 

I 

Tofal 

2.0 

2.0 

2.0 

I 

3.7 

4.2 

14.4 

II 


Table 2: The correlated (type I) and uncorrelated (type II) systematic uncertainties (in percent) on the measure¬ 
ment of the average transverse momentum {pj) of inclusive J/y/ in three centrality intervals in Pb-Pb collisions at 
= 2.76 TeV. 


The non-prompt J/i/r fraction has been determined using an unbinned two-dimensional log-likelihood fit 
described in detail in 071] . which is performed by maximizing the quantity 


N 

InL = In [/s.g -Es.gW •Msig(mee) + (1 -/sig) •MBkg(mee)] , (5) 

where N is the total number of OS candidates in the range 2.2 < mge < 4 GeV/c^ and x is the pseudo¬ 
proper decay length of the candidate 


c(L-pt/pt)'Wj/^ 

x=-(6) 

Pt 

Here L is the vector pointing from the primary vertex to the J/i/r decay vertex and niy^ the mass of the 
J/Y taken from Es,g(.^) and EBkg(.r) (Ms,g(mee) and are Probability Density Functions 

(p.d.f.) describing the pseudo-proper decay length (invariant mass) distribution for signal and back¬ 
ground candidates, respectively. Esig(.x) is defined as 


Es,g(x) = /' • Eb(x) + (1 -/') • F„,{x ), 


(V) 


where Eprompt(.^) and Eb(x) are fhe p.d.f. for prompf and non-prompf J/i/r, respecfively, and /g is fhe 
fraction of reconstructed 1/ i/r coming from beaufy hadron decays 


/b = 




4" ^prompt 


( 8 ) 


A correcfion due fo differenf average (A x e) values, in a given px inferval, for prompf and non-prompf 
1/ i/r, is necessary fo obfain from /g fhe fracfion of produced non-prompf 1/ /b 


/b 



1-/b (Ax£)b 
fs ^ prompt 



(9) 


The various ingredienfs for fhe deferminafion of /b are described in fhe following: 

Monte Carlo pj distributions and polarization assumptions 

Assuming bofh prompf and non-prompf J/i/r to be unpolarized, at a given pj their acceptance times 
efficiency values (A x e) are fhe same. However, fhe pj disfribufions of prompf and non-prompf J/y 
can be differenf, resulfing in differenf average (A x e) compufed over a px range of finife size. Differenf 
hypofheses for fhe kinemafical (px) disfribufions of bofh prompf and non-prompf J/i/r are considered, i.e. 
including or excluding shadowing or suppression effecls as, e.g., fhose predicfed in references 11591-16ill 
for non-prompf J/Y- Due fo fhe weak px dependence of (Axe), fhe resulfing uncerfainfy on /b is small, 
being ~ 5% af low px and ~ 3% in fhe highesf px bin, and is independenf of cenfralify. 
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Fig. 3: The invariant mass (upper panel) and pseudo-proper decay length (lower panel) distributions for e+e pairs 
with px > 1.5 GeV /c in Pb-Pb collisions in the centrality interval 10-50% at = 2.76 TeV. The projections of 
the maximum likelihood fit used to extract /b are superimposed to the data. 


At a given pj, prompt and non-prompt J/i/r can have different polarization and therefore a different 
acceptance. However, the polarization of 1/ i/r from b-hadron decays is expected to be small due to the 
averaging effect caused by the admixture of various exclusive B J/y + X decay channels. Indeed, in 
more elementary colliding systerns, the sizable polarization, which is observed when the polarization axis 
refers to the B-meson direction ll62ll . is strongly smeared when calculated with respect to the direction 
of the daughter J/i/r i63n . as observed by CDF 16411 . The central values of the fraction of non-prompt 


1/ 1 /r are evaluated with Eq. 
as predicted by EVTGEN 152 


assuming unpolarized prompt J/y and a polarization of non-prompt J/y 
. The assumption of a null polarization for non-prompt J/y results in a 
relative decrease of /b by only 1% at high pj (4.5 - 10 GeV/c) and 3% at low px (1-5 - 4.5 GeV/c). The 
relative variations of /b expected in extreme scenarios for the polarization of prompt J/y was studied 
in iHl . The uncertainties related to the polarization of prompt and non-prompt J/ip are are not further 
propagated to the results. 


P.d.f. for prompt J/ Y' Ppmmpt(x) 

The X distribution Epiompt(-^) for prompt J/i/r, which decay at the primary vertex, coincides with the 
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resolution function R{x), which describes the accuracy by which v can be reconstructed. It also enters 
in the p.d.f. describing the v distributions of non-prompt J/i/r, Fij(v), and of the background candidates, 
EflkgW- The determination of R{x) is based on the same MC data sample as used for the inclusive 
J/Y analysis (see Sect. I2.1I) . The systematic uncertainty on R{x) was estimated with a MC approach 
by propagating the maximum observed discrepancies of the track parameters (space and momentum 

and was found to be at most 10%. To 


variables) between data and MC to the x variable II45L1651 
propagate this systematic uncertainty to the final results the fits are repeated after modifying in the log- 
likelihood function the resolution to (l/(l-|-5)) x/?(v/(l-|-5)). In this expression 5 parameterizes the 
relative variation of the RMS of the resolution function and is varied between —10% and -|-10%. The 
systematic uncertainty due to the resolution function is smaller in the highest pj bin, because of the 
better resolution in the x variable and the higher values of the signal-to-background ratio. 

P.d.f. for non-prompt J/Y- Eb(v) 


The shape of the v distribution of non-prompt J/i/r is estimated by using PYTHIA 6.4.21 1671] in the 
Perugia-0 tune ObSH to generate beauty hadrons at ^/s = 2.76 TeV, and the EvtGen package 15211 to 
describe their decays. The systematic uncertainty related to this shape is estimated by assuming a softer 
Pt distribution for the non-prompt J/i/r which is obtained by adding the suppression effects as predicted 
in i 6 in and a harder one taken from the same PYTHIA event generator at ^/s = 7 TeV instead of 2.76 TeV. 
The resulting systematic uncertainty is within 3 - 4%. 

P.d.f. for the background: pBkg (-’C) 

The main difference of the analysis presented in this section, with respect to previous work on pp col¬ 
lisions BTII . concerns the description of FBkg(v). In this analysis such a function includes an extra sym¬ 
metric exponential tail (oc g^bl/'^ym) and depends on the invariant mass and the pt of the dielectron 
pair. It is determined, for each centrality class, by a fit to the data in three pj regions (1.5 - 3, 3 - 4.5, 

4.5 - 10 GeV/c) and in four invariant mass regions on the side-bands of the J/i/r mass peak (2.2 - 2.6, 

2.6 - 2.8, 3.16 - 3.5, 3.5 - 4 GeV/c^, labelled with the indices 1, 2, 3 and 4, respectively), for a total 
of 3 X 4 combinations. The background function in the invariant mass region 2.8 - 3.16 GeV/c^ and in 
each of the three px ranges are obtained by an interpolation procedure as the weighted combination of the 
p.d.f. determined in the other four invariant mass regions. The weights are chosen inversely proportional 
to the absolute difference (or its square) between the mean of the invariant mass distribution in the given 
mass interval and that in the interpolated region 


^Bkgjfiterp kg ,•(''-)> Wi°^ (^ee)interpl (p — 1 Or2). 


( 10 ) 


1=1 


Optionally, only the two adjacent mass regions can be considered in the interpolation procedure, cor¬ 
responding to the condition w\ = W 4 = 0. The central value of /b has been determined as the average 
of the values obtained with the different assumptions (n = 1 or n = 2 , with or without the condition 
wi = W 4 = 0). The RMS of the distributions of the relative variations obtained for /b is used to dehne 
the systematic uncertainty. It becomes larger for central events and in the lowest px interval, where the 
signal-to-background ratio S/B is lower. This approach allows to cope with the much lower S/B ratio in 
Pb-Pb than in pp collisions. 

P.d.f. for the invariant mass distribution of the signal: Ms^g(mee) 

The shape of the invariant mass distribution for the signal is determined by the same MC simulations de¬ 
scribed in Sect. 12. II The influence of defector maferial budgef is sfudied wifh dedicafed MC simulafions, 
where fhe maferial budgef is varied wifhin ifs uncerfainfy (± 6 %) Ibill . The resulfing confribufion to fhe 
sysfemafic uncerfainfy on /g slighfly increases for cenfral evenfs, and ranges from 2 fo 4%. 

P.d.f. for the invariant mass distribution of the background: Ms^g^mee) 

The shape of fhe invarianf mass disfribufion for fhe background candidafes is defermined from ME pairs. 
The relafed sysfemafic uncerfainfy on /g is evaluafed using fhe like-sign disfribufion, insfead of fhe ME 
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one. The uncertainty increases at higher centrality and in the lowest pj interval due to the decrease of 
the S/B ratio. 

As an example in Fig. [3] the projections of the best fit function for n = 1 and w\ = W 4 = 0 are shown 
superimposed to the invariant mass (upper panel) and x (lower panel) distributions of the candidates in 
the centrality range 10-50% for 1.5 < pj < 10 GeV/c. 

A summary of the systematic uncertainties on the determination of the non-prompt 1/ i/r fraction is pro¬ 
vided in Tab.[3]for the three centrality intervals in the integrated pj range and, in the two pj ranges where 
the results will be given, for the most central collisions ( 0 - 10 %). 

The value of /b is determined in two pj bins (1.5 - 4.5 and 4.5 - 10 GeV/c) for the 0-50% centrality 
range and in three centrality classes (0-10%, 10^0% and 40-90%) for 1.5 < pj < 10 GeV/c. The 
/b measurements are then combined with the nuclear modification factors of inclusive J/i/r to get the 
non-prompt and prompt J/y Raa 

/■Pb—Pb 1 fPb—Pb 

pnon—prompt J/i//’ 7 b pincl. J/y/ pprompt J/y/’ ^ 7b pincl. J/y/’ /I 

^AA — Tpp ^AA 5 ^AA — Z Tpp ^AA • ^ 

7b ^ ~ Jb 


pp interpolation 

The value of /b in pp collision at ^/s = 2.76 TeV, /b^, is needed to compute the /?aa for prompt and 
non-prompt 1/i/r mesons, see Eq. (fTTl) . It is determined by an interpolation procedure. Therefore, a fit is 
performed to the existing measurements of as a function of px in mid-rapidity pp collisions at ^/s = 
1 TeV (AEICE |37], ATEAS |35] and CMS If^ i. The function used to fit the data is chosen as 


/-°del(pT) = 


J^FONLL 

dpT 


'do, 


phenom. 


J/¥ 


dpT 


( 12 ) 


which is the ratio of the differential cross section for non-prompt J/i/r obtained by an implementation 
of pQCD calculations at fixed order with next-to leading-log resummation (FONEE) 117IH to that for 
inclusive J/i//, parameterized by the phenomenological function defined in Eq. (I3]l. A similar fit is then 
performed to the CDF results i57n in pp collisions at = 1.96 TeV. Finally, the /b^(pt) value at = 
2.76 TeV is determined by an energy interpolation, which gives /b^ =0.122±0.010 in the integrated px 
range 1.5-10 GeV/c. The quoted uncertainty includes: (i) a component from the fit procedure, which 
depends on the uncertainties of both data and FONEE predictions; (ii) the systematic uncertainty due 
to the energy interpolation, which has been estimated by considering different functional forms of the 
y/i dependency (linear, exponential and power law); (Hi) an additional systematic uncertainty, which has 
been obtained by repeating the whole fitting procedure after excluding, one at a time, the data samples 
used for the /b fit in pp collisions at y /5 = 7 TeV. 


3 Results 


Figure |4] shows the (px) of inclusive J/i/r for the three analyzed centrality intervals. The numerical 
values for (px) are summarized in Tab. |4l As a reference, the (px) in pp collisions at the same centre- 
of-mass energy, as determined by the interpolation method described in Sect. 12.11 is also presented. The 
(px) for Pb-Pb collisions is significantly smaller than that for pp collisions. Such a behaviour is not 
observed at smaller centre-of-mass energies (see left panel of Fig. IDl, for which no significant system 
size dependence of (px) is seen. This might indicate the onset of processes which either deplete the high 
Px region or enhance the J/i/r production at low px in heavy-ion collisions at the EHC. The latter effect 
would be expected as a consequence of a significant contribution from cc coalescence. 


It has been suggested i77h that the observable Taa = (Px)aa/(Px)pp should be particularly sensitive to 
medium modifications affecting the J/i/r transverse momentum distributions. The measured (pj) values 
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Fig. 4; The average transverse momentum (pj) of inclusive ]/\j/ measured at mid-rapidity (|y| < 0.8) in centrality 
selected Pb-Pb collisions (filled circles) and pp collisions (open circles) at = 2.76 TeV as a function of the 
number of participants (Apart)- The uncorrelated systematic uncertainties (type II) are depicted by the open boxes. 
Left panel; A comparison to results obtained by the PHENIX collaboration for Au-Au and Cu-Cu collisions at 
= 0.2 TeV 1^ 72] (open and filled diamonds) and by the NA50 collaboration for Pb-Pb collisions at = 
17.3 GeV S (crosses). The (pj) values are calculated for NASO and PHENIX in the pj interval 0-5 GeV/c, 
while for ALICE the pj interval is 0 - 10 GeV /c. Right panel; (pj) is compared to theory predictions by Zhou et 
al. [74] and Zhao et al. |75, T^l for the pj interval 0-10 GeV /c. 



1.5 <pt < lOGeV/c 

Centr. 0-10% 

Type 

Source 

Centr. 

Centr. 

Centr. 

Pt 

Pt 



0-10% 

10-40% 

40-90% 

1.5-4.5 GeV/c 

4.5 - 10 GeV/c 


Resolution function 

23 

15 

10 

28 

12 

I 

P.d.f. for non-prompt J/lfr 

4 

3 

3 

4 

3 

I 

P.d.f. for the background 

22 

15 

5 

23 

19 

II 

MC Pt distribution 

5 

5 

5 

5 

3 

I 

P.d.f. for the invariant mass of signal 

5 

3 

2 

5 

3 

I 

P.d.f. for the invariant mass of background 

7 

5 

3 

7 

5 

II 

Total 

34 

23 

13 

38 

24 



Table 3: Systematic uncertainties (in percent) on the measurement of the fraction /b of J/y/ from the decay of 
beauty hadrons, for different centrality intervals in the transverse momentum range 1.5 < px < 10 GeV/c, and 
in the two pj intervals for the most central collisions. The conttibutions which are fully correlated between the 
different centrality classes are denoted as type I, the uncorrelated ones as type II. 


for Pb-Pb collisions at = 2.76 TeV are summarized in Tab. |4l The corresponding Taa values as a 
function of (Apart) are shown in Fig. |5]and are found to be significantly below unity. This is in contrast 
to results from lower centre-of-mass energies, where either values consistent with unity (PHENIX at 
= 0.2 TeV lH, 221) or around 1.5 (NASO at = 17.3 GeV |73]) were obtained (see left panel 
of Fig. [5]). The measured (Apart) dependences of (pj) and Taa are compared with a transport model 
for inclusive 1/ y/ by Zhao et al. |l75|, 1761] in the right panels of Figs. |4] and [S] This model includes 
regeneration and dissociation processes, based on in-medium J/yr spectral functions, throughout the 
evolution of a thermally expanding fireball. If also incorporates nuclear shadowing by reducing fhe inpuf 
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Centrality 

(px) (GeV/c) 

(Px) (GeVVc^) 

0 -10% 

10-40% 

40-90% 

2.23 ±0.10 ±0.08 
2.01 ±0.12 ±0.08 
2.02 ±0.19 ±0.29 

5.50 ±0.58 ±0.25 

4.97 ±0.65 ±0.34 
5.15± 1.05 ±1.23 

PP 

2.54 ±0.02 ±0.01 

9.07 ±0.15 ±0.07 


Table 4: The numerical values of {pj) and (p^) calculated in the range 0 < px < 10 /c for the three analyzed 
centrality intervals in Pb-Pb collisions (the first uncertainty is the statistical and the second is the uncorrelated 
systematic (type II), the correlated uncertainty has a value of 2%, see Tab. la. The values for pp collisions obtained 
by the interpolation procedure are given as a reference. 



0 100 200 300 400 0 100 200 300 400 


(N ) (N ) 

part parr 


Fig. 5: The ratio Taa = {p^)aa/ (Pt)pp in the pj interval 0-10 GeV/c for inclusive i/xj/ measured at mid-rapidity 
(|y| < 0.8) in centrality selected Pb-Pb collisions (filled circles) at = 2.76 TeV as a function of the number 
of participants (Apart)- The uncorrelated systematic uncertainties (type II) are depicted by the open boxes, while 
correlated uncertainty (type I) is shown as the filled box at unity. Left panel; A comparison to results obtained by 
the PHENIX collaboration for Au-Au and Cu-Cu collisions at = 0.2 TeV Is, 72] (filled diamonds) and by the 
NASO collaboration for Pb-Pb collisions at = 17.3 GeV if^ (crosses). The PHENIX and NASO Taa values 
are calculated in the pj interval 0 - S GeV /c. Right panel; Taa is compared to theory predictions by Zhou et al. 

1741 and Zhao et al. [75, ^ for the pj interval 0-10 GeV/c. 


charm cross section by a factor of up to 1/3, with a centrality dependence as estimated in 1178h . There is a 
fair agreement between our (px) results and the model calculation, while the Taa is not described by this 
prediction. Our (px) and Taa results are also compared with the calculations by Zhou et al. 11741] . These 
calculations are also based on a transport approach and incorporate dissociation and regeneration of J/i/r 
and heavier charmonia, as well as nuclear shadowing according to EKS98 JH. While the most central 
data point is matched by the prediction, it does not describe the evolution of Taa towards peripheral 
collisions. It must be noted that our results from Pb-Pb collisions at forward rapidity iF/il] exhibit a 
continuous decrease of (px) and Taa from peripheral towards central events and are thus closer to the 
theory predictions, while the behaviour of mid-rapidity Pb-Pb results is more compatible with a flat 
(Apart) dependence. 

The 7 ?aa of inclusive J/i// in three px intervals is shown in Eig.j^along with the results by the CMS col- 
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Fig. 6: The nuclear modification factor /?aa of inclusive J/y/, measured at mid-rapidity (|y| < 0.8) in Pb-Pb 
collisions (0-40% most central) at = 2.76 TeV, as a function of transverse momentum px- The filled symbols 
are placed at the measured px for the given interval. Since for the data point in 4.5 < px < 10 GeV/c (open 
symbol, 0-50% most central) (px) is not available due to the limited statistics, it is plotted at the centre of the 
Px interval. The uncorrelated systematic uncertainties (type II) are depicted by the open boxes, while the correlated 
uncertainties (type I) are shown as the filled boxes at unity. The data are compared to corresponding results by 
PHENIX for Au-Au collisions (0-40% most centtal) at = 0.2 TeV Jot], by CMS for Pb-Pb collisions (0-40% 
most central) at = 2.76 TeV Q, and to predictions by the model of Zhou et al. if^ and Zhao et al 


laboration for the interval 6.5 < px < 30 GeV/c lll3h . both in 0-40% most central Pb-Pb collisions. The 
corresponding numerical values are 0.82±0.11(stat.) ±0.10(syst.) for the interval 0 < px < 2.5 GeV/c 
and 0.58 ±0.06(stat.) ±0.08(syst.) for 2.5 < px < 6 GeV/c, where the systematic uncertainties quoted 
here are the uncorrelated (type II) ones, as listed in Tab. [T] The data point for 4.5 < px < 10 GeV/c 
corresponds to the /?aa value given in Tab. |5] (centrality range 0-50%). Table |5] also contains the R^a 
values for prompt 1/ i/r, which are numerically identical to the ones for inclusive J/ip. The inclusive Raa 
values below px = 6 GeV/c are significantly higher than those measured at higher px, corresponding to 
a decrease of Raa with increasing px, while the high px data point is close to the CMS measurement. 
This Px dependence is similar to the one observed at forward rapidity il2ll . and is in clear contrast to 


the Px dependence measured at lower centre-of-mass energies by the PHENIX collaboration for 


'‘^NN ~ 

0.2 TeV [Qt]. Eigure^also shows the model predictions by Zhou et al. 17411 . The value of the predicted 
Raa is systematically below the measurement and exhibits a px dependence similar to the one in the data. 
The prediction by Zhao et al. 175117611 is close to our result. In both models, the rise of Raa towards 
Px = 0 is due to the dominant contribution from J/ ip regeneration via coalescence. 

The fraction of non-prompt J/ip in the px range 1.5-10 GeV/c is shown as a function of the number 
of participants for the centrality intervals 40-90% ((Apart) = 38), 10-40% ((Apart) = 192), and 0-10% 
((Apart) = 356) in the left panel of Eig. |7] Within uncertainties, no centrality dependence is observed. 
The Px dependence of /u (centrality: 0-50%) is shown in the right panel of Eig. |7] and compared with 
the measurements by CMS in the centrality interval 0-100% and px > 6.5 GeV/c (for the numerical 
values see Tab. |5]l. Our results at low transverse momenta extend the CMS measurements in Pb-Pb 
collisions towards lower px. Also shown are results at mid-rapidity in pp at ^/s = 1 TeV (AEICE 0711 . 
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Fig. 7: The fraction of 1/ly from beauty hadron decays /b at mid-rapidity measured in the pj interval 1.5 < pj < 
10 GeV/c for centrality selected Pb-Pb collisions at = 2.76 TeV (left). The pj dependence of /b at mid¬ 


rapidity for Pb-Pb = 2.76 TeV, |yj/,|/| < 0.8) and pp (-^5 = 7 TeV, |yj 
with measurements by CDF (|yj/y/| < 0.6) |57|, ATLAS (|yj/,|/| < 0.75) 
(right). 


< 0.9) S collisions is compared 
and CMS (bj/^| < 0.9) EUtI 


PT(GeV/c) 

/b(%) 

i?AA(inclusive 1/ t/r) 

i?AA(prompt J/ty) 

f?AA(iion-prompt J/ty) 

0.0- 1.5 

1.5-4.5 

4.5 - 10.0 

10.7±4.8±2.5 

17.0±6.1±2.2 

0.89±0.20±0.21 

0.76±0.09±0.08 

0.38±0.07±0.06 

0.76±0.10±0.08 

0.38±0.07±0.06 

0.73±0.34±0.20 

0.37±0.15±0.09 


Table 5: The numerical values on the fraction of J/iy from beauty hadron decays /b at mid-rapidity and the nuclear 
modification factors /?aa of inclusive, prompt and non-prompt J/ly for Pb-Pb collisions at ^/s = 2.76 TeV. These 
results correspond to the centrality interval 0-50%. The first uncertainty is statistical and the second uncorrelated 
systematic (type II). 


ATLAS ll35h and CMS 17011 ) and in pp collisions at y/s = 1.96 TeV (CDF LSTII ). Considering the ALICE 
and CMS results in Pb-Pb collisions together, a similar pj dependence as in pp is observed. However, 
this similarity could be coincidental, being due to a compensation of the medium effects on the prompt 
component (J/iy dissociation and recombination) and on the non-prompt part (b-quark energy loss). 

In Fig. [8] the nuclear modification factor for non-prompt J/i/r for 1.5 < pj < 4.5 GeV/c and 4.5 < 
Pt < 10 GeV/c is shown together with the result by CMS for 6.5 < pt < 30 GeV/c 113] and with 
theoretical model predictions OOL l3lL 159146iL I80l485ll. One should note that the centrality ranges are 
not the same for AEICE (0-50%) and CMS (0-20% and 20-100%). However, the results obtained by 
CMS for these two centrality bins are compatible with each other, and also compatible with our mea¬ 
surement in the high pj interval (4.5 < pt < 10 GeV/c). The model by Uphoff et al. 116IH follows a 
partonic transport approach based on the Boltzmann equation, which allows interactions among all par- 
tons. It does not include radiative processes for heavy quarks. The calculation has been performed for a 
fixed impact parameter b = 5 fm. In the model of Alberico et al. 1159116011 the propagation of the heavy 
quarks in the medium is described by the relativistic Eangevin equation. The predicted px dependence 
of Raa is strongly influenced by the choice of transport coefficients. Two values are considered, either 
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Fig. 8: The nuclear modification factor /?aa at mid-rapidity (|y| < 0.8) for non-prompt J/ y/ in Pb-Pb collisions 
at = 2.76 TeV as a function of transverse momentum pj. The ALICE measurement corresponds to the 0- 
50% centrality range and to the pj intervals 1.5 < px < 4.5 GeV/c and 4.5 < px < 10 GeV /c. The uncorrelated 
systematic uncertainties (type II) are depicted by the open boxes, while the correlated uncertainties (type I) are 
shown as filled boxes at unity. Results by CMS for higher px in the centrality range 0-20% and 20-100% Q are 
also shown (the two points have been slightly displaced horizontally for better visibility). The data are compared to 
theoretical predictions at mid-rapidity (see text for details). In the right panel, the ALICE result in the px interval 
4.5 < Px < 10 GeV/c is compared to theoretical predictions integrated over the same px range. 


as provided by a perturbative calculation (hard thermal loop approach) or extracted from lattice-QCD 
simulations. The calculations have been provided for the centrality range 0-50%. A transport approach, 
which is based on a strong-coupling scheme, is employed in the model of He et al. 88011 . The trans¬ 
port is implemented using non-perturbative interactions for heavy quarks and mesons through the QGP, 
hadronization and hadronic phases of a nuclear collision. In particular, the elastic heavy-quark scattering 
in the QGP is evaluated within a thermodynamic T-matrix approach, by generating resonances close to 
the critical temperature that can in turn recombine into B mesons, followed by hadronic diffusion using 
effective hadronic scattering amplitudes. The hydrodynamic evolution of the system is quantitatively 
constrained by the measured transverse momentum distributions and elliptic flow of light hadrons. Ra¬ 
diative processes, which should improve the description at high px> are not included in this approach. 
The calculations have been performed in the centrality range 0-50%. The model of Vitev et al. |^, 
assumes the existence of open heavy flavour bound-state solutions in the QGP in the vicinity of the crit¬ 
ical temperature. A description of beauty quark quenching is combined with B meson inelastic breakup 
processes. Furthermore, modified beauty parton distribution functions and beauty fragmentation func¬ 
tions in a co-moving plasma are implemented in this calculation. The prediction is shown for a fixed 
centrality, corresponding to (Apart) = 200, a value very close to the average number of participants in 
the centrality range 0-50%. In the model, a sizable fraction of the suppression is ascribed to the inelas¬ 
tic break-up processes (collisional dissociation), as can be deduced from Fig. [8] by comparing the full 
model prediction with and without the contribution of this specific process. The model of Djordjevic 
ii, shown in Fig. [8]for the centrality range 0-50%, uses a formalism that takes into account finite size 
dynamical QCD medium with finite magnetic mass effects and running coupling. In the WHDG model 
88211 (centrality range 0-50%) the energy loss is computed using perturbative QCD and considering both 
elastic and inelastic partonic collisions and path length fluctuations. The approach of Aichelin et al. 
8831 18411 includes a contribution of radiative gluon emission in the interaction of heavy quarks with light 
quarks, which are considered as dynamical scattering centers. In this model the relative contribution to 
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the energy loss by radiative processes, as compared to collisional ones, is influenced by introducing a 
finite gluon mass. The results of the model shown in Fig. [H which are obtained for the centrality range 
0-50%, correspond to either a pure collisional scenario or a combination of collisional and radiative 
energy loss. Finally, in the model of Horowitz and Gyulassy iSSH . also applied to the centrality interval 
0-50%, the string inspired AdS/CFT gravity-gauge theory correspondence 87] is applied to the case 
of heavy quark energy loss. In the right hand inset of Fig. [H the AFICE 7 ?aa value, integrated over the 
range 4.5 < px < 10 GeV/c, is compared to theoretical predictions computed in the same px range. Most 
of the models predict a larger value of /?aa than observed in the measurement. However, more precise 
data are needed to discriminate among the different models. The next FHC run will provide increased 
statistics for this measurement. 


4 Conclusions 

A study of J/i// production at mid-rapidity in Pb-Pb collisions at = 2.76 TeV has been presented. 
A reduction of the inclusive J/i/r (px) is observed in Pb-Pb collisions in comparison to pp. The ratio 
fAA = (Px)aa/(Pt)pp found to be significantly below unity, corresponding to a medium-induced change 
in the shape of the px spectra. The nuclear modification factor /?aa depends on px. It is around 0.8 for 
Px < 2.5 GeV/c and reaches, at higher px, almost the same level of suppression as observed at RHIC 
energies at low px- These observations might be indicative of a sizable contribution of charm quark 
coalescence to the J/i/r production at low px. Transport models including this additional component are 
able to qualitatively describe the features seen in the data. 

The fraction of J/i/r from beauty hadron decays is determined as a function of centrality and px. No 
significant centrality dependence is observed. By combining this measurement with the inclusive J/i/r 
results the R^a of non-prompt J/yr is obtained in the region 1.5 < px <10 GeV/c, thus extending the 
coverage of CMS to the low px region. The nuclear modification in the region 4.5 < px < 10 GeV/c is 
found to be stronger than predicted by most of the models. 
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